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I. INTRODUCTION 



Branching ratios of B meson two-body nonleptonic decays have been measured by CLEOIII, Belle and Babar 
collaborations [|l]j2[||,^. CP violations in these modes may be observed in near future. Cognizant of this point, 
we have presented some theoretical anticipations for the B Ktt tttt, np and KK decays in the 
perturbative QCD (PQCD) framework. In particular, 5 ~ 15% CP violation is expected in the B Ktt decays. 
The B ^ n, p transition form factors are the integral part of two-body nonleptonic decay amplitudes. In this 
paper we shall convince readers that these form factors in the large recoil region of light mesons are calculable 
in PQCD. This is where our approach starts to differ from other approaches to exclusive B meson decays. 

According to PQCD factorization theorem, a form factor is written as the convolution of a hard amplitude with 
initial-state and final-state hadron distribution amplitudes (f>{x), where x is the momentum fraction associated 
with one of the partons. It has been pointed out that perturbative evaluation of the pion form factor suffers 
nonperturbative enhancement from the end-point region with a momentum fraction a; — > ||^. If this is true, 
the hard amplitude is characterized by a low scale, such that expansion in terms of a large coupling constant 
as is not reliable. More serious end-point (logarithmic) singularities have been observed in the twist-2 (leading- 
twist) contribution to the i? — > tt transition form factor The singularities even become linear at twist 
3 (next-to-leading twist) |ll|]. Because of these singularities, it was claimed that the B ^ n form factor is 
dominated by soft dynamics and not calculable in PQCD We shall argue that this conclusion is false. We 
shall show that at the end points, where the above singularities occur, the double logarithms as In'^ x should be 
resummed in order to justify perturbative expansion. The result, called threshold resummation jT^ , !!^ , leads to 
strong Sudakov suppression at x — + Therefore, the end-point singularities do not exist in a self-consistent 
PQCD analysis. 

In this work we shall investigate contributions to the _B — > tt and B ^ p transition form factors from twist-2 
and from two-parton twist-3 distribution amplitudes. 

In Sec. II we illustrate the PQCD formalism by studying the pion electromagnetic form factor. We review 
the reasoning why one might conclude that the form factor is not calculable, and explain why these objections 
are not justified in QCD. 

In Sees. Ill and IV we derive the B meson transition form factors. It will be shown that the twist-3 con- 
tributions, which seem to be proportional to mo/A/s or Mp/Ms, do not vanish in the Mb — > oo limit. Here 
mo, Mp, and Mb are the chiral symmetry breaking scale, p meson mass, and B meson mass, respectively. We 
record our results of the form factors at large recoil: the B ^ tt form factor F+ ~ 0.3 and the B ^ p form 
factor Aq 0.4. 

Meson distribution amplitudes are defined and the Sudakov factor from threshold resummation is derived in 
the Appendices. 
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II. PQCD APPROACH TO FORM FACTORS 



The suggestion that a hadronic form factor is calculable in PQCD was first made in Ref. [|16|,g/7|,^ . The 
rough idea is summarized as follows. One expands the bound-state wave function for a pion in terms of Fock 
states containing on-shell partons (quarks or gluons) 



*M = i^igq) + i^iqqg) + iPiqqgg) + iPiqqqq) + tpiqqqqg) H — . 

Define a soft function 'J m (A) at a typical hadronic scale A as the initial wave function, 

*m(A) = ^^(qq) + ^P^iqqg) + ^^{qqgg) + i>^{qqqq) + ^^{qmlg) + ' ' ' • 
The wave function \E'm can be related to \1/m(A) via 

*M = *M(A) + G^if*M(A) , 

where K is an irreducible kernel and the Green function involving only hard loop momenta. 
The pion electromagnetic form factor Ft^{Q'^) is then expressed as a convolution integral, 



(1) 



(2) 



(3) 



dxidx2(fki^(fk2±ip^{Pi,xi,ki±^)THiPi,xi, fci^; Pi + q, X2, k2i_)^p^{Pi + q, X2, k2i_) 



(4) 



with Pi being the momentum of the initial-state pion, q large momentum transfer, and 



'P' . Here we have 



written the parton momenta associated with the initial state and final state as k\ — [xiQ /2,kii^,xiQ /i) and 
k2 = (a;2Q/2, k2^, — X2Q/2), respectively, in the notation = (jP , , , p^) , and made explicit the dependence 
of the two-parton wave function ip^{qq) = 'ip^{Pi,xi,ki±) on Pi and ki. The first term in Eq. (Q) contains 
leading contributions, and ellipses represent those from higher Fock states, which are down by powers of l/Q^ 
in the light-cone gauge and by powers of as- The leading diagrams are displayed in Fig. 0. It can be shown that 
the large momentum transfer flows through the hard amplitude Th, and that all nonperturbative dynamics 
goes into wave functions. One can therefore compute Th perturbatively. 
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FIG. 1. Leading-order contribution to Ft,{Q'^). 

However, it was pointed out that the above argument suffers a grave difficulty the diagrams in Fig. 1 
may be infrared divergent, because important contribution to the form factor comes from the region, where the 
exchanged gluons are soft. PQCD is then not applicable. Below we shall examine this difficulty in more details. 
According to Eq. (|^), the first diagram in Fig. |l| gives 
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dz d' 



3 ^ \ 



(27r)3 (27r) 



7r{P2)\d^{y)upmO) 



(5) 



with the color factor Cp = 4/3, the number of colors Nc = 3, and the hard amplitude 
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(6) 



Write (fc2 — fci)^ ~ — a;ia;2Q^ — |fci± — fc2±P- If we ignore |fcij^ — fc2j_p, it has been shown that the integral in 
Eq. (j^) is dominated by contributions from the end-point regions with xi, X2 —> 0. If the pion wave function 
does not vanish at x — > 0, the integral will be even infrared divergent. If we somehow regulate the infrared 
singularity by an appropriate choice of the wave function, the running coupling constant as{xiX2Q'^) evaluated 
at the hard gluon momentum xiX2Q^ is still too large to make sense out of the perturbative expansion. 



A. Feynman's picture of a form factor 

The above end-point singularity corresponds to the picture of the pion form factor Feynman had in mind. 
In the so-called brick-wall frame the initial-state pion with momentum Pi = ((5/2, 0, 0, Q/2) is struck by a 
space-like current of momentum q — (0, 0, 0, —Q), and turns around with momentum P2 = ((9/2, 0, 0, —Q/2) as 
shown in Fig. ^. Feynman pointed out that the major contribution to the form factor comes from the region, 
where one of the partons carries the full pion momentum. The rest of partons, being all wee partons, do not 
know in which direction they are moving. The resulting configuration is essentially identical to the initial-state 
pion except that the momentum of the fast parton is reversed. Hence, Feynman claimed that the dependence 
of the pion form factor is related to the probability of finding a single parton carrying all the pion momentum. 
Feynman's picture is consistent with the statement that the form factor is dominated by the singular part of 
Eq. (H). Because it is singular, we can not compute the form factor. 



Pj= (Q/2A0,Q/2) 



n 



A 



I 



q = (0,OA-Q) 



A 



P2= (Q/2A0,-Q/2) 



FIG. 2. Feynman's viewpoint of the dominant contribution to the pion electromagnetic form factor. 



3 



We argue that Feynman's picture of the pion form factor is false. Consider a QED example. When an 
electron undergoes hard scattering, it can not help but emit infinitely many photons in the direction of the 
electron momentum. As a consequence, the elastic scattering cross section ^(e^e^ — + e^e^) at finite angle 
vanishes at high energy, implying that the probability for the final e^e~ state being accompanied by no photons 
diminishes. In other words, the final state must be accompanied by many photons. In the QCD case of the pion 
form factor, when a quark inside the pion gets hit by a current, the final state will contain many gluons unless 
the spectator quark is nearby to shield the color charge. When one of the quarks carries all the momentum, the 
rest of the pion can not shield the color charge of the fast quark, and many gluons will be emitted in arbitrary 
directions during the hard scattering. Thus, the final configuration ending up as a single pion is extremely 
unlikely. This is so-called Sudakov suppression on exclusive processes at kinematic end points. Therefore, the 
contribution from Fig. ^ is negligible, and the end-point singularity does not exist! 

In the above argument we have ignored the term |fci^ — fc2_Lp in Eq. (|^). Small in momentum space 
corresponds to large transverse distance b of two valence quarks. The color charge of the quark, which is struck 
by the current, is not shielded in this large b region, and will emit many gluons. The probability for having 
a single pion in the final state is then vanishingly small. That is, this configuration can not contribute to the 
form factor. Hence, the momentum space with k± ~f 0, where the end-point singularity occurs, is also Sudakov 
suppressed jl^. The typical behavior of the Sudakov factor exp[— 6, Pi)], x = 1 ~ xi, which is associated 
with the struck quark, is shown in Fig. ||. We observe that the Sudakov factor decreases fast at large b for 
a; '--^ 1 (xi ~ 0), which corresponds precisely to the end-point region in Eq. (|^). In conclusion, the end-point 
singularity is absent, and the major contribution to Fig. || comes from the region with hard gluon exchanges. 




FIG. 3. The Sudakov factor exp[—S{x,b, Pi)]. Note that its value is very small in the region h ~ femax ~ 1/Aqcd, 
with the QCD scale Aqcd = 250 MeV. 



B. Twist-3 contributions 



As derived in the Appendix A, a light-cone pion distribution amplitude is written as 



{7T-{P)\d^iy)upmO) = --= / dxe"P-y {[75 ^^y^A^) + bsl/S^^^o^W 

V Zl\c Jo 

+mo[75(/i+ A- - ^)]i3',€ix)} , (7) 



where P — (P+, 0, 0^) is the pion momentum, the Ught-Uke vector z — (0, z~ ,0±) the coordinate of the d quark, 
and the dimensionless vector n+ = (1,0, Oj^) paraUel to P and n_ = (0, 1,0^) parallel to z. Here a four-vector 
has been expressed in terms of light-cone coordinates, 

The distribution amplitude 0,^ is twist-2, and (ftl^ and (/)^ proportional to toq — M^/ (md + mu) ^ 1-4 GeV, where 
niq is the current quark mass of the quark q, are twist-3. The origin of these terms can be simply understood 
by means of the field-current identity from chiral symmetry, 

d'jBU = imofj.n^ . (9) 

It is easy to observe that twist-3 contributions are suppressed by a power of itiq/Q. The asymptotic behaviors 
of (pTT, 4'lr ^'^d 4>P are known to be 

(j)„{x) (X x{l - x) , 4>P'*ix) (xl . (10) 

As the hard amplitude in Eq. (|^) is convoluted with these distribution amplitudes, we find that twist-2 con- 
tribution is finite, while twist-3 ones are logarithmically divergent without Sudakov suppression. The Sudakov 
factor then introduces an effective cut-off to the integral at Xc ~ Aqcd/Q, and the twist-3 contributions are 
proportional to (toq/Q) ln((5/AQCD)- That is, the power counting is not altered by a logarithmic divergence 
in the factorization formula. As shown later, the power counting for contributions to the B meson transition 
form factors is modified by linear divergences in the factorization formulas. Therefore, the different end-point 
behavior leads to different power counting rules for the pion form factor and for the B meson transition form 
factors. 



III. B TRANSITION FORM FACTORS 



In the B meson rest frame, we define the B meson momentum Pi and the pion momentum P2 in the light-cone 
coordinates: 

Mr Mr 

Pi = ^(1,1,0^), P2 = ^(^,0,0J, (11) 

with the energy fraction 77 carried by the pion. The spectator momenta ki on the B meson side and k2 on the 
pion side are parametrized as 

/ Mb - \ / Mb - \ 

ki = \^0,xi—-^,ki^j , k2 = \^X2ri—^,0,k2j_j . (12) 

Note that the four components of ki should be of the same order, 0(A), with A = Mb — rUb, rrib being the 
b quark mass. However, since k2 is mainly in the plus direction with /c^ ^ 0{Mb), the hard amplitudes will 
not depend on the plus component k^ as explained below. This is the reason we do not show kf in Eq. ([l^ ) 
explicitly. 
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Consider the configuration for the semileptonic decay B irlv depicted in Fig. ^ which corresponds to soft 
contribution to the B tt form factor F^"^ . The u quark and the lepton pair fly back to back with energy of 
Mb /2. The spectator quark d carries a momentum of 0(A). If this configuration is responsible for the decay, it 
is impossible to compute F^"^ using PQCD. However, applying an argument similar to that used for the pion 
form factor, we know that the u quark recoiling against the lepton pair is bound to emit infinitely many gluons. 
Thus, Fig. ^ in fact corresponds to the inclusive decay B X^lv. The probability that the final state in Fig. ^ 
contains only a single pion is suppressed by the Sudakov form factors. A quantitative estimate of Sudakov 
suppression of the soft contribution to F^'^ in the QCD sum rule formahsm will be discussed later. 




FIG. 4. Soft contribution to F^''. 



A. Threshold and k±_ resummations 

It has been explained that the internal b quark involved in the hard amplitude becomes on-shell as the 
momentum fraction x of the d quark vanishes ||l5[| . The contributions to the B ^ tt form factor F^'^ are then 
logarithmically divergent at twist 2 and linearly divergent at twist 3. We argue that as the end-point region is 
important, the corresponding large double logarithms In^ x need to be organized into a jet function St{x) as 
a consequence of threshold resummation p^ . This jet function vanishes as x — s- 0, 1, and modifies the end-point 
behavior of meson distribution amplitudes effectively. This modification provides a plausible explanation for 
the model of the twist-3 pion distribution amplitude proportional to x{l — x), which was adopted in Q. Our 
numerical study shows that the results of the B —i- n form factor obtained in this work are almost the same as 
those obtained in j^. In the following analysis we shall employ the approximate form, 

S,{x) = ^-^^^^^^^^[xil-x)r, (13) 
VttI (1 + c) 

where the parameter c « 0.3 comes from the best fit to the next-to-leading-logarithm threshold resummation 
in moment space. Note that the jet function St is normalized to unity. For details of the derivation, refer to 
the Appendix D. 

Similarly, the inclusion of kj_ regulates the end-point singularities, and large double logarithms as In^ k± are 
produced from higher-order corrections. These double logarithms should be also organized to all orders, leading 
to k± resummation pO| , ^ . The resultant Sudakov form factor, whose explicit expression can be found in our 
previous works ]2^ , p3t , controls the magnitude of k'j^ to be roughly 0{AMb) by suppressing the region with 
k'j^ ~ O(A^). The coupling constant Q;s(AA/B)/7r ^ 0.13 is then small enough to justify the PQCD evaluation 
of heavy-to- light form factors fl. We emphasize that the hard scale for heavy-to-light decays must be AMb 
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in order to define a gauge-invariant B meson distribution amplitude [g^. We shall include the Sudakov factor 
associated with the light spectator quark of the B meson. Whether this factor is essential will be determined by 
the B meson distribution amplitude. Since the B meson is dominated by soft dynamics with xi ^ 0{A/Mb), 
the associated Sudakov effect is minor compared to that from the energetic pion. 

With the possible order of magnitude of fc^ ^ 0{AMb), a Taylor expansion of the hard gluon propagator 
near the end point, 

1 ^ Zl + + (14) 

is certainly not appropriate. A more reasonable treatment is to keep /c^ in the denominators of internal particle 
propagators, and to drop fc^ in the numerators, which are power-suppressed compared to other 0{AIg) terms. 
Under this prescription, the Sudakov factor from kj_ resummation can be introduced into PQCD factorization 
theorem without breaking gauge invariance of the hard amplitudes. For the same reason, the terms proportional 
to fci ^ 0(A) in the numerators should be neglected. It is then obvious from Eq. (|lj) that the hard amplitudes 
are independent of the component . The dependence of the B meson wave function can then be integrated 
out |24j, leading to the parametrization in Eq. (p^). 

Note that the mechanism of threshold and k± resummations is similar with the former responsible for sup- 
pression in the longitudinal direction and the latter for suppression in the transverse direction. As shown below, 
both twist-2 and twist-3 contributions are well-behaved after including threshold and k±_ resummations. Hence, 
the contributions to F^"^ from Fig. ^ dominate in the large recoil region. In this configuration the d quark gains 
a large momentum parallel to the u quark momentum by exchanging a hard gluon with the 6 or u quark. 



w w 
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FIG. 5. Leading-order contribution to F^'^ . 



B. Form factors 



(^1 + ^2)fj, 2 



F,^ie>Ml^,^ , (15) 



We compute the B ^ tt form factors F+ and Fq defined by the following matrix element, 
{7r{P2)\b{0h,M0)\B{Pi)) = F+{q' 
where q — Pi — P2 is the lepton-pair momentum. Another equivalent definition is 

{7r{P2)\b{0h^u{0)\B{Pi)) - /i(g2)Pi^ + hiq^)P2^ , (16) 
in which the form factors /i and /2 are related to and Fq by 



F+ = l{fi + f2) 



The factorization formula for the _B — > tt form factors is written as 

-12 



(17) 
(18) 



{TT{P2)\b{0M0)\B{P,)) ^ g^CpN, [ dxidx2d'h^d^k2^%^%^e-''--y{n{P2)\d^{y)up 



(0)|0) 



xe''=--(0|&,(0)d5(z)|i?(Pi)) T 



7/3;a(5 



(19) 



The pion distribution ampUtude (7r|(i-y(?/)it^(0)|0) has been supphed in Eq. (|^), and the B meson wave function 
is given by (see the Appendix C) 



(2^)3 ' • " ' ^" ^2Nc 

Employing Eqs. (0) and ([20|), we derive, from Eq. (19), 



[(/>! +Afi3)75<?!'B(fcl)]a5 ■ 



(20) 



h = l&TiMlCFr^ j dxidx2 J hdhhdh(l)B{xi,h)[cl)P{x2) - (j)l{x2)]E{t^'y)h{xi,X2,h,b2) , (21) 
/2 = IGttM^Cf J dxidx2 J bidbib2db24'Bixi,bi) 



(t>j,ix2){l + X2V1) + 2r^ ( X2)(t)iix2) - a;205^(x2) 

77 



+2r,0P^(^(2)/^(a:2,Xl,62,6l)| , 
with the ratio — uiq/AIb and the evolution factor 



E{t^'^)h{xi,X2,b,,b2) 



(22) 



(23) 



In the above formulas we have dropped the terms proportional to the momentum fraction xi ^ 0{A/Mb) as 
argued before, which are power-suppressed compared to the leading terms such as 1 + X2/TI in the form factor 
f2- The explicit expressions of the Sudakov exponents 5b and 5^ are referred to The hard function is 

written as 



h{xi,X2,bi,b2) = St{x2)Ko (^/x^x^jMBbi) 

X [9{bi - b2)Ko (V^Ms^i) /o {y/^MBbs) 
+e{b2 - bi)Ko {V^MBb2) lo iV^MBbi)] , 



(24) 



where the factor St suppresses the end-point behaviors of the pion distribution amplitudes, especially of the 
twist-3 ones. The hard scales t are defined as 



t^'> = max(7^AfB, 1/61, 1/62) , 
t(^) = max(7xIr7Mi3, 1/&1, 1/&2) ■ 



(25) 
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It is obvious that if turning off threshold and k± resummations with as fixed, Eqs. ( ^l]) and (22) are infrared 
divergent. 

We argue that the two-parton twist-3 distribution ampHtudes (/)^'* , though proportional to the ratio toq /Mb , 
need to be taken into account. As stated above, the corresponding convolution integrals for the B — > tt form 
factor are linearly divergent without including Sudakov effects. These integrals, regulated in some way with 
an effective cut-off Xc ^ K/Mb, are proportional to the ratio Mb/^- Combining the two ratios mo/ Mb and 
Mb/^, the twist-3 contributions are in fact not down by a power of 1/Mb- 

mp dx2 ^ /mo\ 



and should be included in a complete leading-power analysis. We emphasize that the presence of linear diver- 
gences modifies the power counting rules, making the difference between the B meson transition form factors 
and the pion form factor. 

Various computing methods have been proposed for the evaluation of the B ^ n transition form factors 
F^'^ [q^] in the literature, such as the lattice technique |2^], light-cone QCD sum rules and PQCD 

[^,^. Obviously, lattice calculations become more difficult in the large recoil region of the light meson. 
However, this region is the one where PQCD is reliable, indicating that the PQCD and lattice approaches 
complement each other. This complementation will be explicitly exhibited in Fig. ^ below. In light-cone sum 
rules, dynamics of the _B — > tt form factors have been assumed to be dominated by the large scale of O(rnh). 
This is the reason twist expansion into Fock states in powers of \/mf, applies to the pion bound state. If this 
assumption is valid, PQCD should be also applicable to the B ^ n form factors. Besides, large radiative 
correction to the B meson vertex, which reaches 35% of the full contribution, or about half of the soft (zeroth- 
order) contribution, has been noticed. This 0{as) correction renders the sum rule for JbF^'^, with Jb being the 
B meson decay constant, quite unstable relative to the variation of input parameters ]2^ , |2^ . To stabilize the 
sum rule, one considers another sum rule for Jb at the same time, which also receives large radiative correction 
to the B meson vertex. The two large vertex corrections then cancel in the ratio /bF^'^/Zb- However, the 
radiative correction to /b is then large. 

A careful look at the light-cone-sum-rule analyses indicates that the soft contribution is more sensitive to the 
end-point (x — > 0) behavior of the pion distribution amplitude than the 0{as) correction [2^ ]. Hence, if the 
end-point behavior of the pion distribution amplitude is modified by the Sudakov factor in this work, such that 
the end-point contribution is not important, perturbative contribution can become dominant. The Sudakov 
effect on the soft contribution to F^'^{0) has been investigated in the QCD sum rule formalism (without 
twist expansion for the pion bound state). In this analysis, the soft contribution without Sudakov suppression 
was estimated to be between 0.15 (corresponding to /b ^ 190 MeV) and 0.22 (corresponding to /b ^ 130 
MeV). The soft contribution to /bF^^ obtained in p8|] is consistent with the above range. It was then shown 
that the Sudakov effect decreases the soft contribution by a factor 0.4-0.7, depending on infrared cut-offs for 
loop corrections to the weak decay vertex. Therefore, the soft contribution turns out to be about 0.06-0.15. 
Compared with the lattice results F^^{0) ~ 0.3, it is reasonable to conclude that the soft contribution amounts 
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to about 30%, which is consistent with the observation made in ||2^. It is a fair opinion that the estimate of 
soft contribution is more model-dependent than perturbative one. For example, perturbative contribution is 
less sensitive to the pion distribution amplitude or to other input parameters such as the Borel mass in light- 
cone sum rules [ p8| . In the PQCD approach we calculate the perturbative contribution to which is more 
model-independent, and show that the result can more or less saturate the value predicted by lattice technique. 
For the B meson distribution amplitude, we adopt the model 

2 .2 L " 



4>B{x,b) ^ Nbx'^{1 ~ x)'^ exp 



1 /xMsV 

2 \ LJB J 



(27) 



with the shape parameter lob = 0.4 GeV |5j. The normalization constant Nb is related to the decay constant 
/b = 190 MeV through the relation 

fB 



dxi(l)B{xi,0) 



It is easy to find that Eq. ( |2^ ) has a maximum at x ^ A/ Mb- We employ the models for the pion 

3/. 



07r(a;) 



/2N, 



■.x{l -x) 1 + QAACl'^{2x - 1) + 0.25C4/^(2a; - 1) 



1 + 0A3C^^'^{2x - 1) + 0.09C4'^(2a; - 1) 
(1 - 2x) [1 + 0.55(10a;2 - 10a; + 1)] , 



1/2/ 



2V2A^c 

with the pion decay constant = 130 MeV. The Gegenbauer polynomials are defined by 



(28) 

(29) 
(30) 
(31) 



cTit) = ^(5t2 - 1) , Cl'\t) 



1 



^1/2/ 



1 



(35t4 



15 



(21i4 



SOt^ 4 



3) 



(32) 



whose coefficients correspond to mo — lA GeV. 

We first investigate the relative importance of the twist-2 and twist-3 contributions to F+{q^), and the results 
are listed in Table |. It is observed that the latter are in fact larger than the former, consistent with the argument 
that the twist-3 contributions are not power-suppressed. The light-cone sum rules also give approximately equal 
weights to the twist-2 and higher- twist contributions to Fj^ [Q. We then compare our results of Fj^{q^) and 
Fo{q'^) for = ^ 10 GeV^ with those derived from lattice QCD |^ and from light-cone sum rules in 
Fig. ^, where lattice results have been extrapolated to the small region. Different extrapolation methods 
cause uncertainty only of about 5% [Q. The good agreement among these different approaches at large recoil 
is explicit. The fast rise of the PQCD results at slow recoil indicates that perturbative calculation becomes 
unrehable gradually. The values of F+{0) = Fo(0) = F{0) from PQCD for the parameter lub = 0.40 ± 0.04 GeV 
are listed in Table ||. The resultant range i^+(0) — 0.30 ± 0.04 is in agreement with F{0) ^ 0.3 obtained in 
[^,^. We shall adopt the same range of in the evaluation of the B ^ p transition form factors below. We 
also examine the uncertainty of our predictions from the parametrization of the jet function in Eq. (|l3[). The 



values of F^{q^) vary about 15% for the choices of c = 0.2 and c = 0.4 as shown in Table III. The variation for 
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Fo{q'^) is similar. In a future work we shall incorporate the exact jet function into a convolution integrand in 
moment space. 



(GeV^) 
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10.0 


twist 2 


0.120 


0.128 0.138 


0.148 0.159 


0.172 0.188 0.204 


0.223 


0.243 


0.270 


twist 3 


0.177 


0.193 0.210 


0.230 0.253 


0.279 0.308 0.344 


0.385 


0.432 


0.487 


total 


0.297 


0.321 0.348 


0.378 0.412 


0.451 0.496 0.548 


0.608 


0.675 


0.757 


TABLE I. Contributions to F+(g2) 


from the twist-2 and two-parton twist-3 pion distribution amplitudes. 




iVB (GeV) 


0.36 


0.37 


0.38 0.39 


0.40 0.41 


0.42 


0.43 


U.44 


F(0) 


0.:V15 


0..3.31 


0..321 0..309 


0.297 0.287 


0.277 


0.268 


O.zoy 






TABLE II. 


Values of F+(0) = 


Fo(0) = F{0) for given ujb. 








(GeV^) 


0.0 


1.0 2.0 


3.0 4.0 


5.0 6.0 7.0 


8.0 


9.0 


10.0 


c = 0.2 


0.347 


0.376 0.406 


0.442 0.482 


0.527 0.580 0.639 


0.709 


0.790 


0.886 


c = 0.3 


0.297 


0.321 0.348 


0.378 0.412 


0.451 0.496 0.548 


0.608 


0.675 


0.757 


c = 0.4 


0.260 


0.280 0.303 


0.330 0.359 


0.392 0.432 0.475 


0.527 


0.588 


0.659 



TABLE III. Values of for c = 0.2, 0.3, and 0.4. 
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IV. B ^ p TRANSITION FORM FACTORS 

Consider the semileptonic decay B plv in the fast recoil region of the p meson |3^] . We define the B meson 
momentum Pi as in Eq. the momentum P2 and tlie polarization vectors e of the p meson in light-cone 

coordinates as 

eL = ^^{v\-rl,QA_) , er = (0,0, 1,0) or (0,0, 0,1), (33) 

with the ratio Vp = Mp/Ms and the energy fraction 77 carried by the p meson. We first keep the dependence 
of the kinematic variables in Eq. (|33|), and extract the twist-3 terms proportional to Vp. The parametrization 
of P2 and e is chosen to make this extraction straightforward. 

The B —I- p form factors are defined through the following decompositions of hadronic matrix elements. 
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Mb + Mp 

{p{P2,e*)\mn,umB{P,)) = 2MpA^{q^)'—^q>^ + {Mb + Mp)A,(q^) 



e-q 



Ml - Ml 

1 ' 2 ^ 2 



(34) 



(35) 



To calculate the form factors V , Aq Ai and A2, we adopt the following procedures. First, only the transverse 
polarization vectors tr are involved in Eq. (Q) and associated with the definition of Ai in Eq. (|35|), through 
which we evaluate the form factors V and Ai, respectively. Both the structures associated with Ai and A2 are 
orthogonal to the Icpton pair momentum q. Contracting Eq. ( ^5|) with q^, we have 



(p(P2,e*)|&(0) ^75"(0)|S(Pi)) = 2MpAo(g^)e* - q , 



(36) 



which implies that only the form factor Aq is relevant in two-body nonleptonic decays such as i? — > pT:{K). 
We calculate Aq from Eq. ( ^6|) using the distribution amplitudes associated with a longitudinally polarized p 
meson. 

For the longitudinal polarization vector , the structures of Ai and A2 are in fact proportional to each other: 



e-q 



Mb + Mp 



„ Ml -Ml 

plj, ^ B_ 



{e*-qf{Ml-Ml-q^) 
{MB + Mp)[{t* -qf +q^] 



which can be easily derived via the relation, 

Mf, - Ml - q^ 

P1 + P2- ' 



(e* • qY + q^ 

Contracting Eq. ( |35| ) with e*^ — e* • qq^./q'^, we obtain 



^ M^-Ml + ie*^ 



(e* • q)^ + q^ 



(37) 



(38) 



{p{P2,e*)m 



2P2-q {e*-qy 
Mb + Mp gf2 



q- 



A. 



75w(0)|i?(Pi)) 



{Mb + Mp)2 



1 



(39) 



2P2-q 

from which the form factor A2 can be computed. It turns out that Ai and A2 have a simple relation, since the 
left-hand side of Eq. ( p9[ ) is power-suppressed. 

We derive the leading-power factorization formulas. 



V — SnAIgCp J dxidx2 j bidhxb2db24>B{xi,bi) 



4>p {X2) + rp\\-+X2] 0°(X2) - X2(j)l{x2) 



E{t^^^)h{xi,X2MM) 



+rp [<j)l{x2) + rp{x2)\ E{t^^^)h{x2,Xi,b2,bi) 

Aq = SttMbCf J dxidx2 J bidbib2db24'B{xi,bi) 



(40) 



(1 -I- iqX2)(l)p{x2) +rp[{l- 2x2Wp{x2) ( - - 1 - 2x2 ) (^'p{x2) 



E{t^^^)h{xi,X2,bi,b2) 
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Ai—SnAIgCFi] J dxidx2 j bidhib2db24>B{xi,bi) 



(41) 



4K,{x2) + rp 



+ X2 C(a;2)-X20"(x2) 



E{t'^^^)h(xuX2,bi,b2) 



-rp [(P;ix2) + rpix2)] Eit^^^)hix2,x,,b2,b,) 



A, 



V 



(42) 
(43) 



with the evolution factor E{t) the same as in Eq. (|23|). Taking the fast recoil limit with ?7 — > 1 and assuming 
the asymptotic behavior 0^ — (pp, the above form factors are found to obey the symmetry relations |pl] , p4| , 



V = Ai, A2^Ai- 2r„^o 



(44) 



where the term — 2rpAoj being higher-power, does not appear in Eq. (43). Note that the form factors, treated as 
nonperturbative objects, are not calculated in Instead, the diagrams we have calculated above are regarded 
as perturbative corrections to the relations in Eq. (|4^). 

We adopt the p meson distribution amplitudes given in the Appendix B |3q| , 



(j)p{x) 



^pix) 



3/p 
J p 



■.x{l -x)\l + O.lSCf ^(2a; - 1) 



{3(2a; - 1)^ + 0.3(2a; - lf[b{2x - if - 3] 



2V2]Vc 

+0.21[3- 30(2a;- 1)2 +35(2a;- l)"*]} , 



2J2N, 



:(1 - 2x) [1 + 0.76(10x2 - \Qx + 1)] 



/2Nc 
fp 



x{l-x) 1 + O.2C2' (2a; - 1) 

"^-[1 + {2x - 1)2] + 0.24[3(2a; - 1)^ - 1] 



2x/2]V; I 4 
'0.12[3 - 30(2a; - 1)^ + 35(2a; - 1)*] 



3/p 



4:V2N, 



:(1 - 2x) [1 + 0.93(10x2 - lOx + 1)] 



(45) 

(46) 
(47) 

(48) 

(49) 
(50) 



with the decay constants fp — 200 MeV and fj — 160 MeV. The dependence of the form factors V and 
^0,1,2 with the same B meson distribution amplitude in Eq. ( p7| ) and Alp — 0.77 GeV employed, is displayed in 
Fig. 1^. Our results are consistent with those from light-cone QCD sum rules at small q'^. 
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It is found that the symmetry relation V = A\ in Eq. (^ij) holds very well: A\ is larger than Y only by 
2% in the large recoil region, even after considering the pre-asymptotic forms of 0^ and (\fp in Eqs. ( ^ ) and 
( ^0| ) , respectively. To compare our results with the second symmetry relation, we include next-to- leading power 
terms in Eq. (py), obtaining 



A. = 



l + 2r. 



Ai 



V 



(1 + r]X2)4>pix2) + r 



dxidx2 I t'id6i62'^^20s(2;i, 6i) 
3 



2V 



- 1 cb'M + (1 - 2X2^X2) 



E{t<-^^)h{xuX2,bi,b2) 



+2rp(biix2)E{t<-^^)hix2,x,,b2,b,) 



(51) 



Because of the cancellation of the term 2rpAi/r] and the second term in the above expression, the values of A2 
only slightly deviate from those in Eq. (^). The numerical study shows that A2 is larger than Ai — 2rpAQ by 
about 40%, which can be regarded as the estimate of the symmetry breaking effect. 



V. CONCLUSION 



In this paper we have presented a complete leading-power and leading-order PQCD evaluation of the B ^ tt, 
p transition form factors in the large recoil region. It has been shown that under Sudakov suppression arising 
from fcj^ and threshold resummations, the end-point singularities (logarithmic at twist 2 and linear at twist 3) 
do not exist. The soft contribution to the form factors, being Sudakov suppressed, becomes smaller than the 
perturbative contribution. The physical picture for the mechanism of Sudakov suppression has been discussed. 
We have emphasized that the twist-3 contributions are in fact not power-suppressed in the Mb 00 limit. The 
treatment of the parton transverse momenta k± and the light spectator momentum ki in the B meson in the 
computation of the hard amplitudes has been clearly explained: the hard amplitudes should not be expanded 
in powers of fc^ as the end-point region is important. Using the light meson distribution amplitudes derived 
from QCD sum rules, and choosing an appropriate B meson distribution amplitude, we have derived reasonable 
results for the B tt, p form factors, which are in agreement with those from light-cone QCD sum rules and 
from lattice calculations. Our study indicates that in a self-consistent perturbative analysis, the heavy-to-light 
form factors are calculable. 

The jet function from threshold resummation needs more thorough exploration. We shall investigate the 
relevant subjects, such as factorization theorem in moment space, threshold resummation up to next-to-leading 
logarithms, application to nonleptonic B meson decays and numerical effects elsewhere. Note that if 
considering only /ct resummation twist-3 contributions, though infrared finite, are still too large to give 
reasonable heavy-to-light transition form factors, because the large double logarithms as In^ x have not yet been 
organized. 



16 



Acknowledgements 

We thank S. Brodsky and H.Y. Cheng for helpful discussions. The work was supported in part by Grant-in 
Aid for Special Project Research (Physics of CP Violation), by Grant-in Aid for Scientific Research from the 
Ministry of Education, Science and Culture of Japan. The work of H.N. L. was supported in part by the National 
Science Council of R.O.C. under the Grant No. NSC-89-2112-M-006-033, and by Theory Group of SLAG. The 
work of T.K was supported in part by Grant-in Aid for Scientific Research from the Ministry of Education, 
Science and Culture of Japan under the Grant No. 11640265. 



APPENDIX A: PION DISTRIBUTION AMPLITUDES 

It has been shown [24| that the factorization in fcrmion flow between the pion distribution amplitude and the 
hard amplitude is achieved by inserting the Fierz identity, 

Ii]Iik = 2^iklij + ^(75)^/0(75);^ + j{lf.)ikh^)ij + I (757^)^^(7'' 75)ii + ^{cr^.u)^k{<J'"")lj , (Al) 
into the quark and anti-quark lines of the pion, where / represents the identity matrix, and a^'^ is defined by 



(jt^t^ — i[7^7'^ — j'^^'^]/2. The insertion of Eq. (Al) then leads to various nonlocal matrix elements 



(0|u(0)757prfWk-(P)) , {0\u{Qh5d{z)\7T~{P)) , {0\u{0)^5<J^,d{z)\TT- (P)) , ••• (A2) 

each of which is characterized by different twists. The light-like vector z — (0, z^ , 0±) is the coordinate of the 
d quark, and P = (P+,0, Ox) the pion momentum. 

The general expressions of the relevant matrix elements are, quoted from 

(0|S(0)7M75d(2)k-(P)) =*^P^ f\xe-'-''-'M^) 



(0|S(0)75d(^)|vr-(P)) =-*^mo / dxe-'-"--'' cj,p{x) , (A4) 

^^c Jo 

(0|u(0)75a^.d(z)|7r-(F)) = ^^mo (l - ^\ (P^z, - P,z^) f dxe-"^-^0,(:r) , (A5) 

where (f) and g-^ are the distribution amplitudes of unit normalization, Af^ the pion mass, x the momentum 
fraction associated with the d quark. It is easy to observe that the contribution from independent of the 
pion mass, is twist-2, and the contribution from is twist-4 because of the factor M^. The contributions 
from (pp and (f>a, proportional to = mo/ Mb, arc twist-3. We shall neglect the twist-4 terms and the term 



(A'U/moy in Eq. (A5) 



It is straightforward to read off the pseudo-vector and pseudo-scalar structures of the pion distribution 



amplitudes from Eqs. ([A^) and (A_4). To derive the pseudo-tensor structure from Eq. (Af;), we need more effort. 



Using integration by parts, Eq. (A5) is rewritten as 
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{0\u{0h5<7^..d{z)\TT-iP)) = i^mo (^1 - j\xe-^-P-^^cj,,[x) 



(A6) 



with the anti-symmetric tensor e^j/, = 1. The tensor e^^^ in Eq. (A6) contracts to the spin structure 
a^^^^j2 in the evaluation of the corresponding hard amphtude. The factor 1/2 comes from the extra factor 1/2 
associated with the pseudo-tensor structure compared to other structures in Eq. (|Al[). We have 



T.<^i,v(y^''ih = "2(7^7 - 7 7"^)75 = /»+ - 1)75 



Therefore, up to twist-3, the initial-state tt meson distribution amplitudes are written as 

»i 



(0|«(0),d(z),|7r-(P)) =- 



2A^c Jo 



(A7) 



(A8) 



with 



07r(a;) 



2x/m 



2x/2iV^ 



For the final-state tt meson, we consider the adjoints of Eqs. (A_3), (A4) and (|A5|) 



Jo 



(7r-(p-)|J(z)a^,75w(0)|0) 



6iV, 



Too 1 



2 / "^M" 







ax 



(A9) 

(AlO) 
(All) 
(A12) 



It is observed that the pseudo-tensor structure in Eq. (A12) acquires an extra minus sign, compared to the other 
two structures. The pseudo-tensor structure is then given by — 75(/i- /i+ — 1) = 75 /i- — 1). Therefore, 
up to twist 3, we have Eq. for the final-state tt~ meson. Note that there is an extra term in the definition of 
which contains a differential operator applying to hard amplitudes [0. This term, being power-suppressed, 
is negligible here. The distribution amplitudes (jj^^ and 0^ are normalized according to 



^"'^^^"^ ^ A 



dxdfjx) ■ 
' 2V2iVe 



The tensor distribution amplitude is normalized to zero, because of 



dx — 4)a{x)=4>a{l)-^o{0)^Q 

Q dx 



if vanishes at the end points of the momentum fraction. 



(A13) 



(A14) 



APPENDIX B: p MESON DISTRIBUTION AMPLITUDES 



We choose the p meson momentum P with P^ = M^, which is mainly in the plus direction. The polarization 
vectors e, satisfying P • e = 0, represent one longitudinal polarization vector cl and two transverse polarization 
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vectors ct- Their explicit expressions in light-cone coordinates have been given in Eq. (|33|). To arrive at the 
factorization in fermion flow, we insert the Ficrz identity into the quark and anti-quark lines of the p meson. 
The spin structures in Eq. ( |Al| ) lead to the following nonlocal matrix elements, 

(p-(F,e)|d(z)7^u(0)|0) , (p-(P,e)|J(z)a^,u(0)|0) , 

{p-{P,e)\d{z)Iu{m) . (p-(P,e)M»7p75«(0)|0) , (Bl) 

characterized by different twists. The definition of z is the same as that for the pion distribution amplitudes in 
the previous Appendix. 

The general expressions of the above matrix elements are, quoted from ]35| ], 

-\^,j^r:^Ml j\xe^^-[g,{x)-<j,\\{x)]^ , (B2) 
{p-{P,e)\diz)a^MOm^~i^{{eT^.Pu~eTuP^) I dxe^P-'M^) 



-{P^z, - P,z^)^-—^M'^ I dxe-^-/if|^(x) 
1 . . 



1 M f 

+ -ieTt.z, - ^Tuz^)jr^ dxe'-P-^hsix) - M^)] > , (B3) 



(p-(P,e)|J(z)/.(0)|0) = (^/J - /p^^^^^j e ■ zMlj^ dxe^^^-h\;\x) 

1 /,T , m^ + mA e-z 2 „P.. d (,) 



_ _ /„ ••"" ' — IM,^ / dxe'^''-' — h) \x) , (B4) 
2Nr\ P -"^ Mp J P ■ z P .L dx ^ ' ' ^ ' 



{p-(P,e)\d{z)TolMm) = {j^P-J^p '^^M^' ) Mpe;:-PeT.P^^zp j\xe^^P-g^^\x) , 

- 4^^(^/p /p j p^^e^.pae^P n_ dxe ^^3t W > (ti5) 

where fp and fj are the decay constants of the p meson with longitudinal and transverse polarizations, respec- 
tively, and X the momentum fraction associated with the d quark. We adopt the convention e"^^'^ = 1 for the 
Levi-Civita tensor eA"'a/3_ ^pj-^g distribution amplitudes </>, g and h arc normalized to unity. 

Following the similar procedures, we derive the p meson distribution amplitudes up to twist 3, 

1 



(p-(P,ei)|d(z),u(0),|0) = -== / dxe"^-{Mp[/i]y0p(x) 

\/Zl\c Jo 

+ [^L n.^lix) + Mp[I]i,^;{x)} , (B6) 
{p-{P,eTMz)MOm = I dxe^^^-{Mp[Mi,^;{x) + \h /P\iAl{.x) 



c Jo 

+j^^^,^pAl5l%ei^P''n^^rp{x)'j , (B7) 
for longitudinal polarization and transverse polarization, respectively. We have dropped the terms proportional 



to (twist-4) and the terms (tok -|- nidi/Mp in Eqs. (B4) and (p35|). The definitions of the above distribution 
amplitudes are, 
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(B8) 
(B9) 



APPENDIX C: B MESON DISTRIBUTION AMPLITUDES 



(2^)4 



According to | |39|] , the nonlocal matrix element associated with the B meson is written as 

-^e'''^-'{OK{0)ds{z)\B{P,)) 
{^i+Mb)iz 
{/Pi+Mb)To 



i 



/2Nr 



(j}B{ki) 



V2 



(t>Biki) 



with the wave functions, 



(CI) 



(C2) 



Because the light meson momenta have been chosen in the plus direction, the hard amplitudes for the heavy- 
to-light transition form factors are independent of the component as explained in Sec. III. We construct the 
B meson distribution amplitude, 



(/)(a;i,6) = / dk+d^ki±e'^'^-''(l3{ki) , 



(C3) 



with Xi = k^ /P^. 

The two B meson distribution amplitudes 4'b{x) = (/>^(x, 0) and 4>~jj(x) = (tyjj{x^ 0) are related by the equation 
of motion p9[. 



(C4) 



Assuming that 0^ vanishes at both ends of the momentum fraction, x — s- and x — > 1, we derive 

fB 



1 .1 

dx(f)]^{x) = / dx(j)g{x) 



dxx(j)g{x) — 2 I dxx(j)g{x) 



2^/2Nl 

2A Jb 



(C5) 



/o JO Mb 2y/2Nl 

Therefore, (j)B is normalized to zero. 

We shall argue that the contribution from the distribution amplitude 0b is negligible compared to that from 
(j)B- Consider the reasonable parametrizations, 



(j)B{x) 

4>Bix) 



fE 



2^/2Nc 

fB 

2^f2Wr 



b\x- 



A 



A 

Ml 



2Mb V Mb) \Mj 



2M 



6' X- 



A 

Ml. 



O 



A^ 



(C6) 
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whose moments satisfy Eq. (C5). As shown m Sec. Ill, the hard amplitudes are approximated by ln(l/xi) at 
small xi. A simple estimation indicates that the contribution from (j)B , proportional to ln(MB / A) , is numerically 
larger than that from (pB, proportional to a constant. Hence, after taking into account Eq. ( |C4| ), we consider 
only a single B meson distribution amplitude in this work. 



APPENDIX D: THRESHOLD RESUMMATION 



In this Appendix we supply details of the derivation of the Sudakov factor in Eq. ([T^ ) . Threshold resummation 
introduces a jet function St{x) into the PQCD factorization of the B —>■ it form factors near the end points p^ , 

ra+*oo ^]Sl St{N) 



St{x) 



2-Ki N 



-{l^x) 



-N 



(Dl) 



where a is an arbitrary real constant larger than all the real parts of poles involved in the integrand. The factor 
\/N comes from Mellin transformation of the initial condition sf\x) — 1, 

r-l 



dxil~xf-'si'\x) = ^ 



(D2) 



The Sudakov factor St{N) in the moment (N) space has been derived explicitly to the accuracy of leading 
logarithms (LL) 

1 



- {LL) , 2 i\r 

-7}^ ^In N 



with the anomalous dimension 7]^^'' — olsCfIt^- The contour integral in Eq. (Dl) leads to 



Sf''\x) = -exp (Ja.Cf) F ^(1 - xT"^^"^ sin (^asCj.t') exp (-^C^i^ 



(D3) 



(D4) 



which vanishes at a; 0, since the integrand is an odd function in t, and at a; ^ 1 due to the factor (1 — 2;)°'^^*^*) . 

In this paper we consider threshold resummation up to next-to- leading logarithms. At this level of accuracy, 
the anomalous dimension 7^ contains two-loop contributions, and the coupling constant becomes running. 
The Sudakov factor St{N) is then given by 

1 fi-i/^ dz /-(i-^)^ dX 

2 Jo 



St{N)^exp 



I — Z .l(l-z) ^ 



■7k(«.(AA/|/2)) 



(D5) 



with 



TT V TT / 



Ca 



67 7r2\ 


5 


36 12 J 
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V 



(D6) 



Uf being the number of quark flavors and Ca = 3 a color factor. The anomalous dimension is the same as 
that for kr resummation | |40|] . 

It can be shown that St{x) still vanishes at the end points a; — > and x ^ 1. To simplify the analysis, we 
propose the parametrization, 
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Stix) 



2i+2=r(3/2 + c) 



[xil-x)Y 



V^r{i + c) 

whose end-point behavior satisfies the above requirement. In the 



(c 



effects, Eq. (D7) approaches unity. MeUin transformation of St{x) gives 

Sp\N) _ 21+2^(3/2 + c) 



-B{c+l,c + N) . 



(D7) 

0) hmit, i.e., without QCD 
(D8) 



N ^r{i + c) 

The variable N should be large enough to justify threshold resummation up to the next-to-leading logarithms 
a^lniV, and small enough to avoid the divergent running coupling constant as{Mg/{2N'^)) in Eq. ( |D5| ). Per- 
forming the best fit of Eq. ( p^ ) to St{N)/N for 3 < < 7, we determine the parameter c = 0.3. The difference 
St{N)/N - S{'\N)/N is shown in Fig. | for c = 0.2, 0.3 and 0.4. Equation implies that threshold resum- 
mation modifies the end-point behavior of the meson distribution amplitudes, rendering them vanish faster at 
a; 0. 



0.02 
0.01 



--c=0.4 




N 



-0.02 



FIG. 8. Difference between the jet function and its parametrization in the moment space. 
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